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ABSTRACT 

The transient 500 s X-ray pulsar MAXI J 1409-6 19 was discovered by the slit cameras aboard MAXI on Oc- 
tober 17, 2010, and soon after accurately localized by Swift. We found that the source position was serendip- 
itously observed in 2000 during BeppoSAX observations of the Galactic plane. Two sources are clearly de- 
tected in the MECS: one is consistent with the position of IGR J14043-6148 and the other one with that of 
MAXI J1409-619. We report on the analysis of this archival BeppoSAXIMECS observation integrated with 
newly analyzed observation from ASCA and a set of high-energy observations obtained from the offset fields 
of the BeppoSAX/PDS instrument. For the ON-source observation, the 1.8-100 keV spectrum is fit by an ab- 
sorbed power law with a photon index T = 0.87^jo, corresponding to 2-10 and 15-100 keV unabsorbed fluxes 
of 2.7 x 10~ 12 and 4 x 10 -11 erg cm" 2 s" 1 , respectively, and a 2-10 keV luminosity of 7 x 10 34 erg s" 1 for a 
15 kpc distance. For a PDS offset field observation, performed about one year later and showing a 15-100 keV 
flux of 7 x 10 -11 erg cm" 2 s , we clearly pinpoint three spectral absorption features at 44, 73, and 128 keV, 
resolved both in the spectral fit and in the Crab ratio. We interpret these not harmonically spaced features as 
due to cyclotron resonances. The fundamental energy of 44±3 keV corresponds to a magnetic field strength at 
the neutron star surface of 3.8 x 10 12 (1 +z) G, where z is the gravitational redshift. We discuss the nature of the 
source in the light of its possible counterpart. 

Subject headings: X-rays: binaries — Stars: individual: MAXI J1409-619 — Stars: individual: 
IGR J 14043-6 148 
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1. INTRODUCTION 

On October 17, 2010 lYamaoka et alJ (1201 Obi) reported on 
an outburst from a new source, named MAXI J 1409-6 19, re- 
vealed by the Gas Slit Camera (GSC; iMihara et al] 120021) of 
the M onitor of All-sky X-ray Image (MAXI; Matsuo kaet alJ 
2009) experiment aboard the International Space Station. The 
4-10 keV flux decreased from its initial value of ~40 mCrab 
on October 17 to —30 m Crab on the followin g day. 

On October 20 Swift dGehrels etal.1 12004) began a target 
of opportunity o bservation of the MAX I J 1409-6 19 error cir- 
cle (0.°2 radi us). IKennea et alJ d2010bl) used the X-Ray Tele- 
scope (XRT; iBurrows et al.ll2005l) to pinpoint the position of 
the new source with an estimated 90% confidence level un- 
certainty of l."9 and coordinates RA(J2000): 14 h 08 m 02: s 56, 
Dec(J2000): -61° 59' 00."3. Its spectrum is fit by an ab- 
sorbed power law, with photon index T = -0.5±g \. A follow- 
up XRT observation on November 30 found the source —7 
times brighter than on October 2 0, with an average 0.3-10 
keV flux of 7 x 10" ll) erg cm -2 s" 1 (IKennea et alJ2010al) . Only 
in this latter observation a —500 s pulsation was detected, with 
a 42% sinusoidal rms modulation. The source was detected 
since October 18 also in the 15-5 keV energy band by the 
Swift Burst Alert Telesc ope (BAT: iBarthelmv et all 120051) at 
the level of -30 mCrab (IKennea et al.ll2011l) . 

Neither catalogued radio nor X-ray sources were present 
in the Swift error circle, while a 2MASS IR star, the likely 
IR counterpart of the X-ray transient, lays 2"1 from the XRT 



INAF/Istituto di Astrofisica Spaziale e Fisica Cosmica (IASF) 
Bologna, via Gobetti 101, 40129 Bologna, Italy 

2 Dipartimento di Fisica, Universita degli Studi di Ferrara, via Saragat 
1, 44122 Ferrara, Italy 

3 INAF/Istituto di Astrofisica Spaziale e Fisica Cosmica (IASF) Milano, 
via Bassini 15, 20133 Milano, Italy 



position (IKennea et al . 2010b). 

The source was also observed on October 2 2 by the Pro- 
portio nal Counter Arra y (PCA; iJahoda et alJ 119961) aboard 
R XTE ( Bra dt et al.ll993l) . The X-ray spectrum shows a strong 
6.5 keV Iron line (E.W. 200 eV), and a cont inuum modeled 
by refl ection or partially covering absorption ( Yamaoka et al. 
I2010al) . The 2-10 keV flux was about 10" 10 erg cm" 2 s"\ 
and was strongly variable on time-scales of hundreds of 
seconds. An observation performed on December 4 found 
MAXI J 1409-6 19 at a flux level 6-7 times higher than the 
October observation ( Yamam oto et alJ 12010). and confirmed 
the 500 s pulse period observed by Swift. 

From observations performed on December 2 and 3 by the 
GLAST Burst Monitor (GBM; iMeegan etaD 120071) aboard 
the Fermi satellite, the prese nce of a 500 s double peaked pul- 
sating signal was confirmed dCamero-Arranz eFa l. 2010). 

Following the MAXI J 1409-6 19 localizat ion by Swift, we 
searched for BeppoSAX (Boe lla et al.l[l9 97a) archival obser- 
vations with the transient position within the Narrow Field 
Instruments (NFIs) Field of View (FoV). We found that an 
observation performed in the framework of a Ga lactic plane 
surve y contained the MAXI J 1 409-6 1 9 position ( S gue ra et aT] 
l2010h . 

In this Paper we report on the spectral analysis of this Bep- 
poSAX observation, integrated with high -energy data o btained 
from the offset fields of the PDS (Front era et al.ll 19971) instru- 
ment. We further present the analysis of a newly analyzed 
ASCA observation performed during a Galactic plane survey. 

2. OBSERVATIONS 

2.1. BeppoSAX/PDS Observations 

We searched our BeppoSAX archive for observations cov- 
ering the MAXI J1409-619 position. From our local archive 
we are able to extract not only products from the NFI nominal 
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Figure 1. Region centered on the MAXI J1409-619 position (the small black 
circle) with the PDS pointings: ON pointings (red), POFF pointings (green), 
and MOFF pointings (cyan). The PDS FoV is l.°3 (FWHM). Note the two 
PDS observations (one ON and one MOFF) fully covering the source posi- 
tion, and three (one ON and two POFF) partially covering the position (with 
the term "partially" we intend that the source is beyond the collimator FWHM 
but before its FWZI). The bold dotted line corresponds to the Galactic plane. 

pointings (also available from ASDC - the ASI Scientific data 
CenteQ, but also net spectra from the PDS offset fields. 

Indeed, because of the rocking technique used to derive the 
PDS background, two positions offset by 3.°5 with respect to 
the NFI direction were alternatively observed every 96 s. For 
fields free of sources, we expect that differences between the 
count spectra in the "plus OFF" (thereafter POFF) and "mi- 
nus OFF" (thereafter MOFF) positions will be consistent with 
zero. If, on the other hand, a contaminating source is present 
in one of the two offset fields, then the spectral difference cor- 
responds to a net, background subtracted PDS spectrum of the 
offset field. 

We found three sets of observations containing the 
MAXI J1409-619 position in their FoVs. In Figure Q] we 
show a finding chart centered on the MAXI J1409-619 posi- 
tion and showing both the ON-source and OFF-source point- 
ings (the PDS FoV is l.°3 FWHM). In Tabled we detail the 
whole set of BeppoSAX observations used in our analysis. All 
errors are given at the 90% confidence level for a single pa- 
rameter. 

2.1.1. ON-source observations (OP08317 and OP08386) 

The first ON-source BeppoSAX pointing with 
MAXI J 1409-6 19 in the NFI FoV was performed on 
January 29, 2000 within a program aimed at performing a 
systematic study of part of the Galactic plane (OP08317 
- Galactic Plane Survey Field 02). Another ON-source 
observation (OP08386 - Galactic Plane Survey Field 16), 
was performed a week later. This observation partially 
overlaps the OP083 17 field (see Fig. Q). 

The Galactic plane Field 02 net exposures of the two 
imaging instruments Low-Energy Concentrator Spect rome- 
ter (LECS; 0.1-10 keV; 37' FoV; iParmar et all fl997l) and 
Medium-Ener gy Concentrator Spe ctrometer (MECS; 1.8-10 
keV; 56' FoV; iBoella eTaill 1 997bl) are 6812 and 22060 s, re- 
spectively. The difference is due to the constraint that the 
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Figure 2. 2-10 keV MECS image of a mosaic of the two, partially overlap- 
ping, ON-source BeppoSAX observations OP08317 and OP08386, smoothed 
with a Gaussian filter with a a of 24" . Two sources are clearly detected. The 
one present in the MAXI J1409— 619 error box is listed as #2 in Tableff] The 
two "cuts-out" are due to the removal of the calibration source events. 

LECS can be operated only when in the Earth shadow. Good 
data were selected when the instrument configurations were 
nominal, and with an elevation angle above the Earth limb 

>4°. 

A mosaic of the two, partially overlapping, ON-source 
OP08317 and OP08386 MECS images is shown in Fig- 
ure |2l where two relatively bright sources are clearly present 
(see Table |2|. Source #2 position is consistent with 
MAXI J1409-619, while the source #1 posi tion is consisten t 
with the X-ray source IGR J14043-6148 dBird et al]|20"l(ih . 
The third detection (source #3) corresponds to extended low- 
energy emission observed in the LECS image, and may be 
associ ated with a radio sourc e at approximately the same po- 
sition ( Cohen & Green 2 0011). The unc ertainty in the MECS 
positions is 30" dPerri & C apalbi 20021). 

We extracted LECS and MECS data from a 4' (30 pixel) 
radius circular region centered on the MAXI J1409-619 po- 
sition. Concerning LECS data, there are not enough counts to 
allow a spectral reconstruction. On the other hand, we are 
able to extract a 1.8-10 keV MECS spectrum, that we re- 
binned to a minimum of 30 counts per bin to allow the use 
of x 2 statistics. We used a response matrix appropriate for 
the off-axis position of the source. Because the position of 
MAXI J1409-619 is close to the Galactic plane, we checked 
for possible contamination from the Galactic ridge emission. 
To this end, from the MECS image we extracted background 
spectra from two source-free regions, offset by about the same 
angle as the source (to take into account vignetting), and with 
the same extraction radii. These backgrounds are consistent 
with the mean background used routinely for MECS observa- 
tions. For this reason the mean background was used in the 
spectral analysis. 

Given that the PDS is not an imaging instrument, before 
performing joint MECS/PDS spectral fits it is necessary to 
address the problem of determining whether the high energy 
emission observed by the PDS is due to MAXI J1409-619 
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Table 1 

List of PDS pointings covering the MAXI J1409-619 error box 



Seq 


OPn 


Date 


Duration 


Source Name 


Dist 


RA 


DEC 


Type 


PDS rate 


S/N 








(ks) 




(') 


(J2000) 


(J2000) 




(15-100 keV c/s) 




1 


08317 


29/01/2000 


42 


Gal Plane Sur 2 


19.6 


211.3474 


-61.88437 


ON 


0.247 ± 0.069 


3.6 


2 


10480 


07/01/2001 


220 


Circinus Galaxy 


27.3 


212.8952 


-61.80250 


MOFF 


0.450 ±0.032 


14.1 


3 


08386 


07/02/2000 


52 


Gal Plane Sur 16 


49.5 


210.2841 


-61.83828 


ON 


0.086 ±0.049 


1.7 


4 


01652 


23/02/1997 


90 


Alpha Cen 


53.6 


213.5273 


-62.53019 


POFF 


0.086 ±0.052 


1.6 


5 


08484 


19/02/2000 


47 


Gal Plane Sur 23 


70.0 


213.5936 


-61.09281 


POFF 


0.144 ±0.072 


2.0 



Note. — The "Source Name" refers, in case of offset fields, to the nominal ON pointing. The distance listed in column 6 refers to 
the angular distance of the PDS pointing from the MAXI J1409-619 position. 



Table 2 

Sources detected in the MECS field for OP08317 



# Rate 3 pixel RA. Dec S/N Counterpart 

(2-10 keV) x y (20005 



1 45.5±3.6 373.2 317.1 140431.3 -61 47 10.8 12.6 IGR J 14043-6 148 

2 10.6±1.3 188.2 233.2 140800.4 -61 5824.1 8.15 MAXI J1409-619 

3 3.77±0.79 293.9 250.1 14 06 00.5 -61 56 09.6 4.77 



a In units of 10 3 Counts/s 




Energy (keV) 



Figure 3. 1.8-100 keV joint MECS/PDS count rate spectrum (plus signs) 
and the power law best fit model (histogram) of MAXI T1409-619 in the 
OP08317 field. The fit residuals are shown in the bottom panel. 

or to the INTEGRAL source. Bear in mind that, in comput- 
ing fluxes, it is necessary to take into account the triangu- 
lar angular respons e of the PDS collimators (see Fig. 2 in 
iFrontera et aT]|2007l) . 

First, we analyzed the ON-source OP08386 MECS obser- 
vation, that contains IGR J14043-6148 (detected at >15er) 
but not MAXI J1409-619. For this observation we have only 
a marginal PDS detection (see Table [TJ, therefore it seems 
plausible to associate the high-energy emission observed in 
the OP08317 observation with MAXI J1409-619. To fur- 
ther confirm this association, we verified that the PDS spec- 
trum extrapolated at lower energies is consistent with the 
MAXI J1409-619 MECS spectrum. 

In Figure [3] we show the joint MECS/PDS 1.8-100 keV 
MAXI J 1409-6 19 count rate spectrum, together with its 
power law best fit (normalized % 2 of 0.18 for 15 de- 
grees of freedom - dof). It is characterized by a hard 
spectrum, with power law photon index 0.87^ofq, Nh = 
2.8^22 x 10 22 cm" 2 (consistent with galactic absorption, 



Dicke y & Lockmanl 119901) and unabsorbed 2-10 and 15- 
100 keV fluxes of 2.7 x 10" 12 and 4 x 10" 11 erg cm -2 s~\ 
respectively. 

2. 1 .2. MOFF offset field observation (OP 10480) 

From Figure Q] we can see that we have three sets of PDS 
offset observations covering the MAXI J1409-619 position. 
For two of them (OP01652 and OP08484) we do not have suf- 
ficient statistics to perform a useful analysis (see the observed 
count rates in TableQ]). On the other hand, for the MOFF off- 
set field of the Circinus Galaxy observation, OP 10480, data 
have enough statistics for a thorough spectral analysis. 

First, we checked for the presence of any catalogued 
X-ray sources in the MOFF observation FoV. T he only 
source present is the EGRET ( Hartm an et al.l Il999l) source 
3EG J1410-6147, that was rece ntly associated with the 50 ms 
radio pulsar PSR J1410-6132 dO'Brien et alj|2008l) . No sig- 
nificant X-ray emission has been detected around the EGRET 
source (Dohertv et al. 2003); therefore we associate the PDS 
detection with MAXI J1409-619. 

We started by extracting a 15-100 keV coarser binned PDS 
spectrum of the MOFF field (assuming as its background the 
corresponding POFF spectrum). The fit with a power law 
gives a photon index of 2.2 ± 0.3, with a normalized x 2 of 
1.31 for 1 1 dof. From the shape of the residuals and from the 
quite different photon index with respect to the ON-source 
observation, we suspected the presence of a change of slope 
in the spectrum between 20 and 30 keV. This break is not de- 
tected in the ON-source observation because of its poor statis- 
tics. We tried both a broken and a cutoff power law. By fix- 
ing the cutoff energy at 25 keV (we were only able to con- 
strain the break/cutoff energy between 20 and 40 keV), a cut- 
off power law fit gives a power law index T = 0.9±q 3, with 
a normalized \ 2 V °f 1 06 for 1 1 dof. While from a statisti- 
cal point of view the two fits are equivalent (an F-test gives a 
36% probability of chance improvement (PCI) of the \ 2 - see 
Appendix lAl for the discussion on the correctness of the use 
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Figure 4. 15-200 keV MOFF PDS count rate spectrum (plus signs) and the 
best fit power law plus 3 Gaussians in absorption model (histogram). The fit 
residuals are shown in the bottom panel. The Gaussian energies are left free 
and are 44 ±3, 73^, and 128+| keV. The line widths were constrained to 4, 
4, and 7 keV for the three lines. 

of the F-test), the cutoff fit nicely matches the spectral index 
observed at lower energies. The corresponding 15-100 keV 
flux, corrected for the triangular response of the collimator, is 
7.0 x 10" 11 erg cm" 2 s" 1 . 

Intrigued by the shape of the residuals, we tried a smaller 
rebinning factor and a wider (15-200 keV) energy range. We 
were surprised to find the spectrum shown in Figure [4] Ab- 
sorption features at ^40, ^70, and ^120 keV are quite ev- 
ident. We also noted that the points that trail these features 
at ^50, ^90, and ~150 keV seems to be in agreement with 
a hard continuum up to 200 keV. To test this hypothesis, we 
performed a fit to the data below 30 keV together with these 
trailing points. No cutoff (the break energy is not constrained) 
or broken power law (the two photon indices are equal) model 
are able to fit the data, while a power law with photon index 
1 .0 ± 0.2 does (xl = 0.56 for 20 dof). Therefore the change of 
slope in 20^-0 keV derived from the analysis of the coarser 
binned energy spectrum is really due to the presence of the 
absorption features. The first absorption feature drives the 
position of the break, while the higher energy ones are spread 
by the rebinning, resulting in low-count, high energy bins that 
are fit by a change of slope. 

To further support this interpretation, we performed a joint 
fit with the MECS spectrum (this analysis is detailed in Ap- 
pendix IBJ. We find that a cutoff or a broken power law mod- 
els are not able to fit this broad band spectrum with the con- 
straint r ~ 0.9 below 10 keV. This implies that the origin of 
the change of slope is not due to the presence of a cutoff but to 
something else, likely absorption features. For this reasons as 
the continuum we will use a power law with V fixed at 0.87, 
as derived by the analysis on the broad band spectrum of the 
first observation. 

As throughly discussed in Appendix lAl the assessment of 
the statistical significance of these absorption features has to 
be done with great care. Our strategy was to evaluate the 
significance of each feature one at a time. We started with 
a power law fit in the energy range 15-55 keV, and we ob- 
tained an unacceptable xl of 1.4 for 25 dof. The inclusion 
of a Gaussian in absorption (gabs in XSPEC) with centroid 
energy E cyc = 44 ± 3 keV, t = 16+' 4 , and width fixed at 4 keV 
significantly improved the fit (xl equal to 0.48 for 23 dof). 

Two comments on this measurement are in order: first, the 



^44 keV line is statistically significant (PCI is 6 x 10" 3 ). Sec- 
ond, the difficulties in determining the other line parameters 
(the line width had to be fixed, while the line normalization 
is affected by large errors) are due to the fact that the lines 
are seen as "holes" in the continuum: we have not enough 
statistics to clearly reconstruct the line profile. The line width 
is therefore computed from the "hole" width, while for the 
line depth we can infer only an upper limit (this will be more 
evident for the lines at higher energies). 

Following our detection strategy, we performed a second 
fit on the 15-80 keV energy range. The fit with only one 
Gaussian in absorption was not acceptable: xl was 3-6 for 32 
dof. By adding another Gaussian in absorption we obtained 
a significantly better fit, with xl of 0.48 for 30 dof (PCI is 
2 x 10~ 7 ). The second line parameters are E cyc = 73^ keV, 
T = 91^42, while the the line width was fixed at 7 keV. Note 
how the two line energies are not in an harmonic ratio. 

Then we extended the energy range to the whole 15- 
200 keV band. Again the fit with two lines was not acceptable 
(xl of 4.5 for 46 dof), so we added another absorption line, 
the width of which was fixed at 7 keV. The fit significantly 
improved (xl of 0.58 for 44 dof, PCI of 10" 10 ), with the best 
fit E cyc = 128^| keV. r was not constrained (this was not a sur- 
prise because, as discussed above, we have no signal in the 
"hole", and therefore we are not able to put a limit to the line 
normalization). 

Finally, we constrained the line energies to be in an har- 
monic ratio. In this case we obtain a xl of 1.22 for 46 
dof. The best fit parameters of the fundamental were E cyc = 
41 ± 1 keV, t = 1 1+7. All the line widths were fixed, as in the 
previous case with the line energies left free, but we were not 
able to constrain the (1:2) and (1:3) line normalizations. 

In order to verify that the observed features were not of 
instrumental origin, and at the same time to better charac- 
terize the cyclotron resonance features (ther eafter CRFs), we 
performed a normalized Crab ratio analysis ( OrlandinH l2004l) 
on the MOFF spectrum, as we successfully employed for 
the detection of CRF s in numerous X-ray pulsars (see, e.g. 
lOrlandini et alj[l998l) . As it is evident from the top panel of 
Figure El the ratio between the MAXI J 1409-6 19 and the 
Crab count rate spectra shows a "hole" in the 38^-4 keV 
range, due to the fundamental CRF To enhance the feature, 
we multiplied the MAXI J1409-619/Crab ratio by the func- 
tional form of the Crab spectrum, i.e. a simple power law with 
photon index equal to 2.1. The result is shown in the second 
panel of Figure|5] Finally, in the lower panel, we show the ra- 
tio between the previous function and the MAXI J 1409-6 19 
best fit continuum power law model, together with a Gaussian 
fit to the fundamental CRF. The Gaussian width was fixed at 
4 keV, while the Gaussian centroid energy is 43+ 2 keV. Please 
note that this value is a lower limit, because it does not take 
into account the intrinsic energy resolution of the PDS instru- 
ment. 

In order to give a quantitative evaluation of the CRF in the 
normalized Crab Ratio we performed a run test (see detailed 
discussion in Appendix |AV Up to ^34 keV the residuals are 
consistent with random fluctuations (/V+ = 12; AL = 9; N r = 13; 
consistent with random fluctuation at the 84%). On the other 

5 Thereafter we will indicate the fundamental resonance as (1:1), meaning 
with this the ordinal number of the resonance with respect to the fundamental. 
With this notation the first harmonic will be indicated as (1:2), the third (1:3), 
and so on. 



BeppoSAX observations of MAXI J1409-619: discovery of CRFs 



5 



io- 2 

lO- 3 

10- 4 
io-« 

io-= 



10-' 
2 
1.5 
1 

0.5 


-0.5 



. 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


i i i | i i . 

i 1 ■ 

- J i i i i f] 




Wffi 




: , , . 1 , , . . 1 , . . , 1 . , , , 1 , , : 



20 30 40 50 

Energy (keV) 

Figure 5. Top: Ratio between MAXI J 1409-6 19 MOFF and Crab count rate 
spectra. Middle: the MAXI J1409-619/Crab ratio multiplied by E" 2 - 1 , the 
functional form of the Crab spectrum. The dotted line corresponds to the best 
fit power law continuum. Bottom: ratio between the former expression and 
the best fit power law continuum, together with a Gaussian fit to the ~43 keV 
CRF. The statistical assessment of the CRF has been evaluated with a run test, 
finding that we can reject the line being random at the 99% confidence level 
(see text and Appendix [A]i°r details). 

hand, in the 34-50 keV band we see a clear structure in the 
residuals. In this case from a run test with N+ =2, N~= 14, 
N r = 2we we can reject the null hypothesis of randomness at 
the 99% confidence level. 

2.2. ASCA observation 

The region around MAXI J1409-619 was observed by 
ASCA on March 2, 1998 for a net exposure time of about 
18 ks. The source is clearly detected, and its spectrum is con- 
firmed to be quite hard, with a best fit power law photon in- 
dex T = 0.1^:1 and an unabsorbed 1-10 keV flux of 3 x 10" 12 
erg cm" 2 s" 1 . The absorption in the direction of the source is 
not well constrained by the fit, with an upper limit of 3 x 10 22 
cm" 2 , consistent with the Swift and BeppoSAX/MECS results. 

3. A VERY RED(DENED) COUNTERPART FOR MAXI 11409-619 



As stressed bv lKennea et alJ (|2010bl). a NIR sour ce belong- 
ing to the 2MASS catalogue dSkrutskie et al.ll2006l) is found 
within 2"1 of the MAXI J 1409-6 19 position as determined 
by the XRT According to that catalogue, this source (labeled 
as 2MASS J14080271-6159020) has NIR magnitudes J = 
15.874±0.086, H = 13.620±0.022 and K = 12.560±0.021. 
No catalogued optical counterpart is present at the position of 
this NIR source, and inspection of the DSS-ITRed digitized 
archival plates^ does n ot show any eviden t optical object at 
that location. Following Mon et et alJ (120031) we can thus place 
a conservative upper limit, R>20, to the /?-band magnitude of 
the optical counterpart of this source. The above magnitudes 
thus indicate that this is an extremely red object. 

Indeed, assumin g the Milky Way extinction law 
(ICardelli et al.1 [1989), we find that the NIR color indices 

6 |http : / / archive ■ eso ■ org/dss/dss] 



of this sou rce are consis tent with those of a late O/early 
B-type star dWegnerl l l994l) with a reddening of Ay » 20 mag. 
This, using the formula of iPredehl & Schmittl d 19951) . implies 
a column density of N# ~ 3.6 x 10 22 cm" 2 , which is consistent 
with that inferred from the Swift and BeppoSAX spectral 
analysis results. This large extinction is also supported by the 
non-detection of the optical counterpart of the source. All 
this points to a likely HMXB nature of this source, hosting 
a heavily absorbed early-type star, similar to several cases 
of INTEGRAL hard X-ray sources identifie d as HMXBs 
at op tical and/or NIR wavelengths (see, e.g., Mase tti et al.1 
l20loh . 

Assuming then Ay«20 mag along the MAXI J1409-619 
line of sight and a B0 spectral type for the companion star 
in this system, we can infer its distance, depending on the 
luminosity class (main sequence, giant or supergiant) of the 
star. U sing t he tabulated absolute magnitudes for this type 
of star dLandl 19921) we find, for these three cases, respective 
distances of ^4.6, ^7.9 and ~14.5 kpc. Given that the large 
absorption found along the line of sight of this source at both 
NIR and X-ray wavelengths is consistent with the Galactic 
one, we consider most likely that the correct distance is the 
largest one. Thus, MAXI J1409-619 is quite likely located 
in the far side of the Galaxy, i.e. in the farthest parts of the 
Sagittarius-Carina arm, and its mass donor star is likely an 
early-type supergiant. 

4. DISCUSSION 

We found five sets of observations containing the posi- 
tion of MAXI J 1409-6 19 in our BeppoSAX/PDS archive, per- 
formed in 1997, 2000, and 2001, and one in the ASCA archive 
(March 1998). In all our observations the source was in a 
low state, with 15-100 keV fluxes in the range ~2-8 mCrab, 
and no spectral variability during the observations. For com- 
parison, an integrated exposure (over 5 years) of 2.4 Ms by 
INTEGRAL/IBIS provides 2a upper limits on the persistent 
quiescent emission of 0.2 and 0.4 mCr ab in the 20—40 an d 
40-100 keV energy bands, respectively (ISgueraet al.ll2010t) . 
When assuming the source fluxe s in outburst as measured by 
Swift dKennea et al.l2010bl.l201 11) and RXTE (lYamamoto et al.1 
2010), from our low state measurement we can infer a dy- 
namic range of 400 in 15-50 keV, and of 300 in 2-10 keV. 

The discovery of CRFs in the low state spectrum of 
MAXI J1409-619, together with its ~500 s pulsations, un- 
ambiguously identify the source as an accreting X-ray bi- 
nary pulsar. Only few sources show multiple CRFs, and 
only two show resonances above the (1:2). Three CRFs 
were observed durin g the 2004-2005 outburst of the X-ray 
pulsa r V0332+53 dCoburn et al.1 |200H iKrevkenbohm et all 
120051: Tsygank ov et al.1 120061) . and five (possible six) CRFs 
were discovered in th e spectrum of 4U 0115+63 during 
its 1999 giant outburst ([S antangel o et al.|[l999l: iHeindl et all 
1999; Ferrigno et al. 2009). In both cases, deviations from 
a pure harmonic ratio among the CRFs were observed, and 
were explained in terms of departures from a classical dipo- 
lar structure of th e magnetic field in the line-forming region 
dNishimurall2005Ll2008l) . 

The magnetic field strength at the neutron star surface 
corresponding to £ C vc = 44 keV is 3.8 x 10 I2 (l+z) G 
dCanuto & Venturall 19771) . where z, the gravitational redshift, 
for a typical neutron star of mass 1.4 M Q and radius 10 Km, 
is about 0.3. When left as free parameters in the fit, we 
found that the best fit CRF line energies do not follow the 
harmonic relation E„ = nE\. A slightly non-harmonicity is 
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Figure 6. Harmonicity of the MAXI J1409-619 CRF line energies. The 
dashed black line corresponds to a linear fit (that is, E n = nE\) to the data, 
while the colored strips take into account relativistic effects, as detailed by 
Eq. {T}- Each strip corresponds to a fixed value of the magnetic field strength 
B12 = B X 10 12 G, and to a range 0.0001-1 for sin 2 0, where 9 is the angle 
between the photon direction and that of the magnetic field. While E\ and 
£3 are consistent with B12 ~ 3.8, £2 is consistent with a magnetic field about 
20% lower. 

exp ected when relati vistic effects are taken into account (see, 
e.g.. lMeszaroslll992l) 



, l+2»(B/B crit )sin / (9-l 1 

E„ = m e c 7T - — — (1) 

sin 9 1+Z 

where m e is the electron rest mass, c the speed of light, 9 the 
angle between the photon and the magnetic field direction, 
and B c rit = 4.414 x 10 13 G is the critical magnetic field strength 
where the cyclotron energy equals the electron rest mass. 

The observed (1:2) and (1:3) ratios of the line energies with 
respect to the fundamental, 1.7±0.2 and 2.9±0.3, cannot be 
explained in terms of Eq. (Q~|i, as it is evident from Figure [6] 
A fit to the harmonic relation, shown as the dashed line, gives 
Ei = 41 ± 3 keV (in agreement with the best fit value E cyc = 
41 ± 1 keV found when imposing the harmonic relation in 
the spectral fit), but with poor significance = 2.27 for 2 
dof). In the same figure we also show the harmonic relation 
from Eq. (fl]i for different values of the magnetic field and the 
angle 9. Taking into account relativistic effects, we found that 
the magnetic field responsible for the (1:1) and (1:3) CRFs is 
about 20% higher than that responsible for the (1:2) CRF. 

Two points are worth noticing: it is always the (1:2) reso- 
nance that shows the larger disagreement with the harmonic 
relation, and this could be due to the f act that this CRF is 
due to pure absorption (Nishimura 2003), while for the other 
resonances other effects, like multi ple scattering and ph oton 
spawning, enter into play (see, e.g JSchonherr et al.l2007L and 
references therein). Unfortunately, because we are not able to 
reconstruct the CRF profiles, we cannot extract more infor- 
mation, like the electron temperature and the geometry of the 
emitting region. 




Figure 7. Error box of AGL J 1 4 1 0— 6143 (big red circle) superimposed on 
the 20-100 keV INTEGRAL/IBIS deep mosaic image (~2.3 Ms exposure). 
The BeppoSAXfMECS position of MAXI T1409-619 is marked by the green 
square. The two small yellow ellipses represent the Fermi 7-ray sources 
2FGL T1409.9-6129 and 2FGL J1413.4-6204, respectively. White contours 
(from 50% to 99%) refer to the EGRET source 3EG T1410-6147. 

Second, at variance with the V0332+53 and 4U 0115+63 
observations, both performed during giant outbursts, our Bep- 
poSAX observations were performed while MAXI J1409-619 
was in a low state. This did not allow the study of the de- 
pendence of the CRF parameters as a function of luminos- 
ity, an important tool for study of the physical conditions in 
the line-forming region ( Mihara et al. 2004; Nakaii ma et alj 
2006: iKlochkov et al.11201 1» ~ 

The likely early-type optical counterpart and the 500 s pul- 
sation makes MAXI J1409-619 a HMXB pulsar. Accord- 
ing to the nature of the secondary star we have two possibil- 
ities: the source is a supergiant fast X-ray transient (SFXT; 
iSgueraet all 120051: iNegueruela etail 12006]) or a Be/HMXB. 
In favor of the former interpretation is the highly reddened 
supergiant as possible counterpart, the typical outburst X- 
ray luminosity of 2 x 10 37 erg s -1 (assuming a distance of 
14.5 kpc), and the dynamic range of more than two orders 
of magnitude (—300) that is typical of the so-called "interme- 
diate" SFXT dSguera et alJI2007l: IClark et al.ll2010h . A gainst 
we have that the source active phase, about two months long 
dUeno et al. 2010; K ennea et alj|201 1[). is signifi cantly longer 
that that typical of SFXT (see, e.g. JSidoli|[2009h . If this were 
the case, then our magnetic field measurement wo uld rule out 
the magnetar nature for SFXTs dBozzo et al.ll2008l) . The ob- 
served properties of MAXI J 1409-6 19 are also in agreement 
with those observed in other Be/HMXBs, like 1 A 1 118-615, 
a 400 s X-ray pulsar with a hard X-ray spectrum (r ~ 1), a 
CRF at —55 keV, and long (tens of years) periods of quies- 
cence interrupted by giant (Type-II) outbursts lasting weeks 
to months, in which the X-ray luminosit y increases by a fac- 
tor -200 (see, e.g. jRutledge et al.ll2007T) . This would put the 
source a factor —2-3 closer. 

Finally, we note that MAXI J1409-619 is located within 
the 0.°5 error box of the unidentified transient MeV source 
AGL J1410— 6143 (see Fig. [7) dis covered on Februar y 21, 
2008 by the 7-ray satellite AGILE dTavani et alj[2009al) dur- 
ing a bright MeV flare lasti ng only about one day dPittori et alj 
120081: lOrlandini etail 120081) . The AGILE large position un- 
certainty makes very difficult the identification of its lower 
energy counterpart responsible for the 7-ray emission. De- 
spite this drawback, it is intriguing to note that the flaring 
source MAXI J 1409-6 19 is the only catalogued hard X-ray 
source above 20 keV (20-100 keV) to be located inside the 
AGILE error box, according to all the available catalogs in the 
HEASARC database. This spatial correlation is further sup- 
ported by a similar transient nature for both MAXI J 1 409-6 1 9 
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and AGL J1410-6143. We also point out that this is not a 
unique case. To date, a few HMXBs have been unambigu- 
ously detected as fl aring MeV sources lasting only a few 
days (ISabatini et al.l 120101: iTavani et all l2009bt lAbdo et al.l 
2009). In addition, there are several other HMXBs proposed 
as best candidate counterparts of unident ified transien t MeV 
sources located on the Galactic plane (ISguera et al.l 120091 
2011; Sguera 2009). For the sake of completeness, we note 
that within the large AGILE error box there are a number of 
high energy MeV sources (see Fig.|7): i) 2FGL J1409.9-6129 
and 2FGL J 14 13. 4-6204 have been reported in the second 
Fermi source catalog (Abdo et al. 2011) as firmly identified 
7-ray pulsars and this unambiguously excludes their associa- 
tion with the transient AGL J1410-6143, ii) 3EG J1410-6147 
is still unidentified although it has been likely associated with 
the 7-ray pulsar 2FGL J 1409 9-6129 dO'Brien et alj|2008h . 
However, Wall ace et al.l ((2000) reported a possible MeV flare 
from 3EG J1410-6147 lasting a few days on November 1991. 
This behavior is at variance with the proposed association 
with the 7-ray pulsar 2FGL J1409. 9-6129 while it is more 
compatible with the flaring nature of AGL J1410— 6143. Fur- 
ther multi-wavelength studies (radio, NIR, X-ray and 7-ray) 
of the sky region are strongly needed and encouraged to shed 
more light on the nature of such high energy emitters. 
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APPENDIX 

ON THE STATISTICAL SIGNIFICANCE OF ABSORPTION FEATURES 



The analysis of the significance of spectral absorption fea- 
tures, and in particular CRFs, has historically been fraught 
with serious misunderstanding, that we attribute to the appli- 
cation of techniques made for emission features to absorption 
ones. This has a very profound impact on the statistical anal- 
ysis, as detailed in the following. 

First, when searching for spectral features, and especially 
absorption features, the binning of data is crucial: indeed a 
narr ow feature can be lost if the data binning is too high (see, 
e.g. JEadie etal1[l97ll page 259), as we showed in the data 
analysis detailed in Section 12.1.21 But before searching for 
the best binning factor, it has become customary to rebin data 
in order to achieve a minimum number of counts in each bin. 
The reason of this procedure is that we must avoid bins with 
few counts if we want to use the % 2 statistics to test for the 
goodness of our fit (in other words, errors must be normally 
distributed). A mini mum 20 co unts per bin filter is sufficient 
to achieve this goal (Cash1 IT979h . 

The application of this criterion to background subtracted 
spectra is correct for instruments that are not background 
dominated (that is, they have an intrinsically low number of 
counts per bin), like a detector on the focal plane of X-ray 
optics. On the other hand, when the net source spectrum is 
obtained by subtracting two high counts spectra this correc- 
tion must not be applied (instruments like BeppoSAX/PDS or 
RXTE/HEXTE are background dominated, with a very high 
number of counts per energy channel). Although the net 
source counts can be very low, they are still Gaussian dis- 
tributed, because resulting from the difference of two nor- 
mally distributed counts. 

Second, if features are present in the source spectrum, they 
(clearly) show up in the fit residuals, and the standard good- 
ness of fit estimator, the minimization of the x 2 statistics 
(Lam pton et al.|[l976h . should (clearly) indicate that the fit is 
not "good" and an additional component is necessary. The in- 
clusion of a new component to the continuum should result in 
a reduction of the reduced x 2 - In order to assess whether the 
improvement of the x 2 is due to chance or it is because the 
new com ponent is signi ficant it is customary to use the F-test 
(see, e.g.. lBarlow|[l989l page 160). 

A very important point has to be raised here, because we 
have two completely different approaches whether we are 
dealing with emission or absorption features. In the former 
case we are dealing, using a Xspec terminology, with an ad- 
ditive component, and we test th e null hypothesis that the co- 
efficient of the new term is zero (Bevington 1969, page 200). 
This is the F-test routine that is incorporated in the recent ver- 
sions of Xspec. 

On the other hand, when dealing with an absorption feature, 
the component is not additive but multiplicative (see routines 
gabs or cyclabs). It is therefore obvious that we cannot 
use the F-test for an additive component, because we cannot 
vary to zero the new component coefficient without modify- 
ing the parameters of the original model. The two compo- 
nents (the multiplicative model and the original continuum) 
are strongly coupled, and variations in one will affect also the 
ot her. Therefore we must use a different F-test, as described 
in lPress et al.l (120071 page 730), that tests the null hypothesis 
that the observed variances comes from the same sample (and 
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Figure 8. Example on how the goodness of fit estimation performed by the 
X 2 statistics not always is able to find the best fit model to the data. In this 
case the linear fit to the 12 data points yields a normalized x 2 of 1, but the 
deviation with respect to the straight line is evident. A run test shows that 
these deviations have only 1.3% probability to be due to random fluctuations 
(see text for details). 

this routine is not incorporated in Xspec). This same test 
must be used to assess if two different models are sta tistical ly 
equivalent or not (and indeed we used it in Section 12.1.21 to 
compare the power law and the cutoff power law fits). It is 
worth noting that the application of the F-test to the statistical 
assessment of CRFs is not affected by the concer ns raised by 
its use with emission lines (Protassov et al. 2002). 

Third, CRFs are very shallow features, with their count 
rates often at the instrument sensitivity level. We therefore ex- 
pect errors to be dominant, and consequently the x 2 could not 
be the best estimator of the CRF significance. We need a test 
that takes into account the structure of the line (an information 
which is lost in the x 2 , because we compute the square of the 
difference between the data and the model). In other words, 
we need to compute what is the probability that a particular 
structure visible in the fit residuals occurs by chance. 

The run test (also kn own as Wald-Wolfowitz test; Barlow 
[l989tlEadie e"taI][T97Tb works on the signs of the deviations, 
that is on the form of the residuals. To better clarify how the 
run tes t works give a look at Figure [8] adapted from iBarlowl 
(1989): when fitting the 12 data points with a straight line 
the normalized x 2 is exactly 1 (likely due to error overestima- 
tion). But it is evident by eye that the fit is not good (indeed 
the data come from a parabolic model). The reason is that if 
the fit were good we should expect that the number of points 
"above" the fitting line should not group together, but should 
be intermixed with points "below" the fitting line (and this 
should be more true as the number of data points increases). 
If, on the other hand, we observe only small groups of data 
with the same "sign" (called runs), this means that our data 
are not randomly distributed with respect to the fitting model, 
but there is an underlying trend. 

As we can see from Figure [8] we have 12 points, 6 points 
"above" the fitting line (let us call them A^+) and 6 points "be- 
low" the fitting line (N-). The number of runs N r is only 3, 
suspiciously small. Indeed the probability of obtaining N r <3 
is 1.3%, telling us that the structure observed in the residuals 
is not due to random fluctuations but to a wrong modelization 
(the linear fit). 
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From this example, and from the real-case analysis per- 
formed on the residuals shown in Figures [5] and [9] should be 
clear that the goodness of fit in the case of data structures in 
the residuals should not be addressed only with the \ 2 estima- 
tor, but it has to be supported by the run test. 

In conclusion, this is a sort of cookbook for a correct eval- 
uation of the statistical significance of absorption (CRFs) fea- 
tures in the spectra of X-ray sources: 

1. Be careful with the binning: although CRFs are usu- 
ally broad features, a too small binning can hide struc- 
tures. When dealing with background dominated data 
(like, for example, BeppoSAX/PDS or RXTE/HEXTE 
net spectra) never rebin data in order to have a mini- 
mum number of counts per bin, otherwise we risk to 

THE CONTINUUM SPECTE 

The analysis performed on the broad band spectrum of the 
pointed MAXI J1409-619 observation yields as best fit con- 
tinuum a power law with index ^0.9 (see section [2.1.1b . On 
the other hand, the analysis performed on a coarser binned 
spectrum of the offset observation indicates the presen ce of 
a change of slope around 2CMI-0 keV (see section 12.1.21 ). Be- 
cause a finer rebinning of this second observation revealed the 
presence of features that can be explained as CRSFs, it is im- 
portant to understand whether this change of slope is due to a 
cutoff in the continuum or is a combined effect of the binning 
and the presence of the absorption features. 

From the pointed observation we know that data below 
10 keV must be described by a power law with T ~ 0.9. To 
take into account this constraint we performed a joint fit of the 
high energy PDS spectrum from the second observation with 
the low energy MECS spectrum from the first observation. 

Although the two spectra come from different observations, 
X-ray pulsars are known to not show any spectral variability 
besides that at low energy due to reprocessing of X-rays by 
the circumstellar material (see, for example, the case of Vela 
X-l or GX301-2). Because our source does not show any in- 
trinsic absorption or emission lines, signatures of reprocess- 
ing, we are confident that our MECS spectrum did not change 
between the two observations. 

First we tried a fit with a power law, modified at low energy 
by photoelectric absorption. All the parameters are left free 
but the Nh (fixed at the 2.8 x 10 22 cm" 2 value obtained from 
the ON-source observation), and the result is shown in panel 
a. of Figure[9] The fit parameters are listed in Table[3] 

The fit is not statistically acceptable, and the slope is driven 
by the low energy part of the spectrum. To fit also the PDS 
data we need something that is able to change the slope at 
about 40 keV. 

Now let us test whether we can explain this change of slope 
in terms of a cutoff in the continuum. We first tried a cutoff 
power law (see panel b. in Figure[9]i: the fit is statistically ac- 
ceptable but the power law index does not match with that of 



lose the absorption feature; 

2. In order to assess the significance of a CRF the F-test 
is perfectly suitable, but we must use the correct rou- 
tine: the F-test routine included in XSPEC is applicable 
only to an additive component (in particular, an emis- 
sion line), while both gabs and cyclabs are multi- 
plicative; 

3. The evaluation of the statistical significance of a CRF 
must be supported by other tests, especially in the pres- 
ence of structures in the fit residuals: the run test is able 
to discriminate whether the structure is due to random 
fluctuations or not. 



MODEL OF MAXI J1409-619 

the low energy part, below 10 keV. The systematic deviations 
from the model (information lost by the \ 2 test, as discussed 
in AppendixlAt can be quantified by means of the run test on 
the MECS residuals (N+ = 8, 7V_ = 6, N r = 3). The probability 
of obtaining < 3 is 0.5%, therefore the structure observed 
in the residuals is not due to random fluctuations but to the 
wrong modelization of the continuum. The same result is ob- 
tained with the broken power law fit (see panel c. in Figure^: 
there is no match to the MECS data. 



Table 3 

Best fit parameters for the joint MECS/PDS spectral fits 





Power Law 


Cutoff 


Broken PL 


Power Law 

+ 3 gabs 


r 

^cutoff 


0.87 ±0.07 


-0 75 m3 
14±4 


0.06^| 
11% 


1 0+O- 3 
1 '"-0.2 


r 2 

Xl (dof) 


2.48 (66) 


0.91 (65) 


z -"-0.6 
0.78 (64) 


0.54 (59) 



These results demonstrate that the change of slope is not 
in the continuum but is due to something else: we made the 
hypothesis that its origin are absorption features. Therefore 
the true continuum model is a power law together with CRSF 
features (see panel d. in Figure^. Because we are not able to 
disentangle the continuum and the CRSFs, due to the broad- 
ness of the features and their low statistics, we need to fix the 
power law index to the value obtained by the broad band fit. 

Just as a test, we performed all the analysis on the CRSFs by 
fixing the power law index at 1.0, 1.1, and 1.2, corresponding 
to the 1.0 ±0.2 value obtained from the fit to the PDS data 
below 30 keV and the points trailing the CRSFs. The CRSF 
parameters are all consistent with each other, demonstrating 
that they do not depend on the particular value of the power 
law index. 
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Energy (keV) Energy (keV) 

Figure 9. 1 .8-200 keV joint MECS (from the pointed observation) / PDS (from the offset observation) spectrum, a. Fit with a power law continuum, b. Fit with 
a cutoff continuum, c. Fit with a broken power law continuum, d. Fit with a power law continuum and three CRSFs. The only model able to fit simultaneously 
the low (below 10 keV) and high energy data is the power law plus three CRSFs. The best fit parameters are listed in Table[3] 



